Experimental study of z resolution in Acousto-Optical Coherence 
Tomography using random phase jumps on ultrasound and light 
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Acousto-Optical Coherence Tomography (AOCT) is a variant of Acousto Optic Imaging (also called Ultrasound 
modulated Optical Tomography) that makes possible to get resolution along the ultrasound propagation axis 
z. We present here new AOCT experimental results, and we study how the z resolution depends on time 
step between phase jumps T^, or on the correlation length Az. By working at low resolution, we perform a 
quantitative comparison of the z measurements with the theoretical Point Spread Function (PSF). We present 
also images recorded with different z resolution, and we qualitatively show how the image quality varies with 
T(^, or Az. © 2013 Optical Society of America 

OCIS codes: 170.3660, 110.7050, 110.7170, 160.5320, 170.3880 



1. Introduction 
Citation: 

Max Lesaffre, Salma Farahi, Franois Ramaz and 
Michel Gross "Experimental study of z resolution in 
Acousto-Optical Coherence Tomography using random 
phase jumps on ultrasound and light" Applied Optics: 
52,5, pp. 949 - 957 (2013) 

Acousto-optic imaging (AOI) [1-5] combines light and 
ultrasounds to measure local optical properties through 
thick and highly scattering media. 

The principle of this imaging method is the follow- 
ing. Ultrasound (US) is applied in the region of interest, 
within the thick scattering sample, making the scatterers 
vibrate. A CW laser illuminates the sample. The vibra- 
tion of the scatterers at the acoustic US frequency ujus 
(2MHz typically) modulates the phase of the photons 
that are scattered by the sample. The light exiting the 
sample contains thus different frequency components. 
The main component (carrier) is centered at the laser 
illumination frequency ujl- It is related to the scattered 
photons, that do not interact with the US. The sideband 
components at frequencies uj±i are shifted from the US 
frequency so that uj±i = ujl ^^us- The sideband pho- 
tons, which result from the interaction between light and 
US, are called "tagged photons" , their weight depends on 
the optical absorption in the region of interest, where the 
US beam is focused. Acousto-optic imaging detects selec- 
tively the tagged photons. The signal is recorded while 
scanning the US within the volume of the sample, and 
thus one gets informations of the local optical absorption 
and/or scattering in 2D or 3D. 

The first experiments were performed with fast single 
detectors to record the modulation of the optical signal 
at the US frequency [1-3,5,6] , but, since the phase of 



the modulation is different for each grain of speckle, the 
detector can only process one grain of speckle. In or- 
der to increase the optical etendue of detection, Leveque 
et al. [4] have developed a camera coherent detection 
technique that processes many speckles in parallel. This 
technique has been pulled to the photon shot noise limit 
by Gross et al. [7] using a holographic heterodyne tech- 
nique [8] able to detect photons with optimal sensitiv- 
ity [9,10]. But, because of the tissue inner motions, the 
frequency bandwidth of the photons that travel through 
living tissue is enlarged (1.5 kHz typically through 4 cm 
of living breast tissue [11]), while the bandwidth of the 
camera coherent detection (1/T) is limited by the ex- 
posure time T. In order to increase this bandwidth, one 
has to decrease T, and thus the detected signal remains 
almost the same. This means that the detected signal 
is essentially independent of camera exposure time T, if 
this time T is longer than correlation time of the sam- 
ple speckle. This point has been experimentally demon- 
strated with turbid liquid [12]. An alternative to match 
the detection and the signal bandwidth while keeping a 
large optical etendue is to use a detection involving a 
detection bandwidth comparable with the signal band- 
width while keeping a large optical 
etendue, is to use a detection schemes involving photore- 
fractive (PR) crystals [13-17]. 

Since the US attenuation is low in tissues, tagged pho- 
tons are generated along the US propagation z axis at 
nearly constant rate. This means that in a continuous 
regime of the US, the AO techniques give nearly no in- 
formation on the location of the embedded objects along 
the z axis. In order to get z information, acoustic pulses 
can be used whether with single detectors [18,19], cam- 
eras [20] and PR crystals [13-16,21,22]. Nevertheless, 
reaching a millimetric resolution with US pulses requires 



1 



a typical duty cycle of 1%, corresponding to the ratio of 
the exploration length within the sample ('^10 cm) and 
the desired resolution 1 mm). This is problematic 
regarding the very small quantity of light that emerges 
from a clinical sample, since a weak duty cycle yields a 
low signal and a poor SNR. 

In a recent publication, Lesaffre et al. [23] overcome 
the duty cycle problem, and get z resolution with a vari- 
ant of the acousto optic technique, called Acousto Opti- 
cal Coherence Tomography (AOCT). In this technique, 
the same random phase modulation is applied on both 
light and US, but the light modulation is applied at re- 
tarded time t + z/cus^ where z is the coordinate of the 
selected zone, and cus the US velocity. The Lesaffre et 
al. experiment has been made with PR detection of the 
tagged photon signal, and, in that case, the exact shape 
of the z instrumental response has been calculated [24]. 

In this paper, we will present new experimental results 
in order to verify this theory by studying how the z reso- 
lution depends on the characteristic time of the random 
phase jumps, T^. We will recall briefly the theoretical 
results of reference [24], and introduce the AOCT reso- 
lution Az. By working at low resolution, we will perform 
a quantitative comparison of the measured z shape, with 
the theoretical one. We will then present images recorded 
with different resolution Az^ and qualitatively show how 
the images quality varies with Az. 

2. The experimental setup 

The experimental setup is presented on Fig.l. It is close 
to the one described elsewhere [23]. An acoustic trans- 
ducer (Panametrics A395S-SU, focal length of 75 mm) 
generates a US beam, which is focused on the sample 
(frequency uus = 2.3 MHz, with a maximum US pres- 
sure Pus — 1.6 MPa at the focal point). The master laser 
(wavelength A = 1064nm, 5 W after amplification) is 
spatially single mode (transverse and longitudinal) and 
vertically polarized. It is splitted in two beams, whose 
amplitude, phase and frequency are controlled by two 
acousto optic modulators (AOMl and A0M2) driven at 
^AOMi — ^AOM2 = 80 MHz. The first beam illuminates 
the sample (area s ^ 1 cm^). Note that the output sig- 
nal beam is totally depolarized because of the multiple 
diffusion events in the scattering sample. The second 
(reference) pumps the photorefractive (PR) crystal. The 
relative power of the reference versus illumination beam 
is adjusted by a half wave plate (A/2) followed by a po- 
larizing beam splitter (BS). The PR crystal is a bulk 
GaAs crystal (1.4 x 1.6 x 1.4 cm^ [23]), oriented into 
an anisotropic diffraction configuration [22,25]. The sig- 
nal beam, which is depolarized, and the reference beam, 
which is vertically polarized, enter onto orthogonal faces, 
respectively (llV^) and (ll\/2), in order to have a grat- 
ing vector along (001). In such a configuration, the ef- 
fective electro optic coefficient is rg// = r4i. A polar- 
izer oriented at 7r/4 relative vertical is positioned after 
the sample, while a horizontally oriented polarizer is po- 
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Fig. 1. AOCT experimental setup. Laser: Nd: YAG laser 
(1 W @ 1064nm); amplifier: optical amplifier 5W-doped 
Yb; 10: optical isolator Faraday; A/2: half- wave plate; 
AOMl, 2: acousto optic modulators; M: mirror; BS: po- 
larizing beam splitter; X : Radio Frequency (RF) ~ 80 
MHz Double Balanced Mixer; A: RF amplifier; GaAs: 
photorefractive GaAs crystal; PLI, 2: linear polarizers; 
Si: Silicon photodiode (0.3 cm^); L in: Lock-In ampli- 
fier; PC: personal computer; PZT: Piezoelectric acoustic 
transductor (~ 2 MHz). H(t): 0,7r random phase mod- 
ulation at frequency uOmod = 3kHz, and with duty cycle 
r = 0.24; ifjuSi V^p- random sequence of phase. 



sitioned after the crystal. This configuration minimizes 
scattering of the intense reference beam by the PR crys- 
tal, which corresponds to the main parasitic light contri- 
bution. For a typical reference beam of power Pr = 650 
mW over an area Sr = 0.32 cm^, the photorefractive 
response time is tpr = 3.4 ms. 

In order to perform the z axis selection of the tagged 
photons, the phases of illumination {^jjpit)) and US 
{i^us{t)) beams are pseudo-randomized. Every time in- 
terval T<j>, we apply pseudo random phase jumps of 
or TT on both illumination and US (i.e. we multiply the 
corresponding sinusoidal signal by either +1 or -1 every 
T^).By using two DBM (Double Balanced Mixers ZAD- 
IH: Mini- Circuits Lab. Inc.), the sine wave signals of 
frequency ujaomi and ujus that drive the acousto optic 
modulator AOMl and the PZT transductor are mixed 
with two correlated rectangular pseudo random phase 
signals e^'^^^^) = ±1 and e^'^^^^^) = ±1. These pseudo 
random signals are generated by two 80 MHz arbitrary 
waveform generators (Agilent 33250A) with a common 
10 MHz reference clock. Because of the finite memory 
{N words) of the generators, the phase signals, which 
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Fig. 2. Plot of correlation function ^^(r) (a), and its 
square ^^(r) (b). The horizontal axis units is either t /T^ 
(for time correlation), or z/Az with = Cf/s- T«i> (for z 
resolution). 



make random jumps at time intervals of T«^>, are peri- 
odic (pseudo random) with periodicity T = NT^^ where 
N = 2^^ typically. In order to characterize the random 
phase sequence, one can introduce the phase correlation 
function [24]: 

g^ir) = (ei^^(^)e-^-^^(^--)) (1) 

_ (^^j'^Usit)^-ji^Usit-T)^ 

Since N 1, the correlation function g^{r) has a 
triangular shape whose Full Width at Half Maximum 
(FWHM) is in time units, and Az = cusT^ in length 
units (see Fig. 2 (a) ). 

In order to choose the z selected zone in the sample, 
the phase ipp is delayed with respect to ipus in a way 
that 

^p{t)=^us{t-T) (2) 

with r = z/cus- 

The signal detection procedure is the same as in ref- 
erence [23]. In order to get a time varying signal on 
the photodiode PD, we apply on the US beam an extra 
H{t) = ±1 rectangular modulation signal at frequency 
^modl"^^ = ^/Tmod = 3 kHz, with duty cycle r. This 
frequency has been chosen to satisfy the condition of va- 
lidity of theory (see Fig. 2 of reference [24]), 

< Tmod < TpR. (3) 

Other modulation frequencies are possible if Eq. 3 con- 
dition remains fulfilled. We have chosen to use here ail 
rectangular modulation, because this modulation is (i) 
easy to make, (ii) yields a quite good Lock-in signal, and 
(iii) has been considered in theory [23,26]. 
Within the interval t = to Tmod we have thus: 

H{t) = +1 for 0<t/Tmod<r (4) 
= -1 for r<t/Tmod<l 

The US beam is H{t) = ±1 modulated by using a second 
DBM in an AM modulator configuration. A US signal at 
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Fig. 3. Photodiode Lock-in acousto optic signal as a func- 
tion of the longitudinal coordinate z for for Az = 10.0, 
6.0, 4.0, 2.0, 1.6, 1.2, 0.8, 0.6, 0.4 and 0.2 cm. 



- 2 MHz is feeded within the LO port of the DBM, while 
the H{t) rectangular signal is feeded within the IF port. 
The H{t) signal is generated by a computer driven sig- 
nal generator, which controls the frequency ujrnodj'^^ and 
the duty cycle r. The ±1 modulated US output is then 
obtained on the LO port. Since the modulation is fast 
with respect to the recording of the tagged photon grat- 
ing on the PR crystal {oOmod ^ ^/^pr)^ the modulated 
component of the index grating that is grooved in the 
PR crystal is proportional to [(1 — 2r)rsinc(7rr)] where 
r is the extra modulation duty cycle [26]. In order to 
maximize this component in the Lock-in detection, the 
duty cycle is fixed at r = 0.24. 

The photodetector (PD) is a silicon photodiode 
(Hamamatsu S2386-8k, area Spo = 0.3 cm^) connected 
to a home-made lOMl] transimpedance amplifier. The 
modulated component of the Photodiode signal at Umod 
is measured by a Lock-in amplifier (EG&G Inc. 7210; 
integration time r c:^ 100 ms). We have stacked all 
the (sample+crystal+photodetector) together in order 
to optimize the detection flux without using lenses. 

3. Quantitative study of the AOCT Point 
Spread Function (PSF) along z 

Lesaffre et al. [24] have calculated the Point Spread Func- 
tion (PSF) of the AOCT signal, and show that PSF(z) is 
simply equal to g'^ expressed in length units. This func- 
tion is plotted on Fig. 2 (b), and corresponds to what 
expected in experiments. In order to measure PSF(2;), 
and to compare it with theory, we have plotted the pho- 
todiode Lock-In Signal SpD as a function of the delay r 
for phases sequences with different correlation times . 
By this way, one can vary the location z = cus^ of the 
selected zone, and the correlation length Az = cusT^- 
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Fig. 4. Experimental profiles (black circles) made with 
Az = 10.0 (a), and 6.0 cm (b). Theoretical shape gl^ 
(red curves) made by adjusting Az in order to get the 
best fit with experiment. Fit yields Az = 10.27 (a) and 
6.09 cm (b). 
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Fig. 5. Photodiode Lock-In signal SpD as a function of 
the PZT location along the x axis. Experimental data 
(cross), and exp(— x^/it;^) Fit with Wg = 10.2 mm (red 
dashed line). 



The experiment is made with an Agar + Intralipid 10% 
phantom sample without absorbing inclusion. The re- 
duced scattering coefficient is fi'^ = 10 cm~^. The sample 
thickness is L^y = 3 cm along the illumination direction y. 
Other dimension are L^xLz = 6. 5x6. 5 cm^. The optical 
power is = 1 W/1 cm^ on the sample, while the ref- 
erence beam on the GaAs crystal is Pr = 600 mW/0.32 
cm^. Acoustic pressure at maximum is Pmax = 1-6 MPa, 
and lock-in integration time is Tc = 100 ms. 

Although the experiments presented is this paper 
use the same setup, they have been done at different 
time, with slightly different laser power and optical ad- 
justments. The parameters P^, Pr, Pmax and Tc may 
thus change from one experiment to the next. Similarly, 
the sample parameters: thickness Ly^ lateral dimension 
Lx^Lz and scattering coefficient jii'^ may have changed, 
for good or bad reasons. 

Figure 3 shows the Photodiode Lock-In signal as a 
function of z for Az varying from 10 cm down to 0.2 
cm. When Az is small (0.6, 0.4 and 0.2 cm), the curve 
maximum is almost proportional to Az, while the width 
remains nearly constant. The shape of the curves cor- 
responds then to the shape, along of the cloud of 
photons that travel through the sample and that reach 
the Photodiode. In the region crossed by the US beam, 
the thickness of the cloud along x and z is about 1.5 
cm i.e. 0.5 X the thickness of the sample. When Az is 
larger than the width of the cloud, i.e. for Az = 6 and 
4 cm, the curve maximum saturates, and the width in- 
creases linearly with Az. In that last case, one can ne- 
glect the width of the cloud, and the Spd{z) corresponds 
to PSF(z). 

To make a first quantitative test of theory [24], we 
have plotted both the measured acousto optic signal 
and the theoretical point spread function (PSF) gi'^. 
We have considered situation where the PSF width 
{Az = lOandQcm) is larger than the thickness of the 
cloud (about 1.5 cm). In such case, the experimental pro- 
file must be close to theoretical PSF gi'^. We has thus 
plotted gi'^ (red curves of Fig. 4) by adjusting Az in or- 
der to get the best fit with the experimental data . The 
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Fig. 6. Colored curves: experimental profiles Spd{z); 
Red curves: theoretical profiles Agi'^{x) exp(— x^/it;^) 
where (g) is the convolution operator. Az is 10, 4, 2, 1, 
0.5 and 0.3 cm. 



fit yields Az = 10.27an(i6.09 (for the gi'^ curves) in ex- 
cellent agreement with Az = 10and6cm (experiments). 

In order to compare the Spd{z) experimental curves 
with theory for any Az, one has to calculate gi'^{z) 
Cz{z), where is the convolution operator, and Cz{z) 
the shape of the cloud of photons along z. Since the 
cloud is symmetric with respect to the photon propaga- 
tion direction the shapes Cx and Cz along x and z 
are the same. We have measured Cx by moving the PZT 
transducer along x with a step motor, and recorded the 
photodiode signal Spd{x) as a function of the PZT lo- 
cation X. Spd{x) is plotted on Fig. 5. 

We have assumed that Cz{x) = Cx{x) = Spd{x)^ 
and we approximate Cx{x) by this gaussian fit Cx{x) 

'^9 with Wg = 10.2 mm (see Fig. 5). Replacing Cz{z) 
by the fit, we have calculated the theoretical profiles 
Agi'^{x) (g) ex.p{—x'^/Wg)^ where A is an arbitrary gain 
parameter that should be the same for all value of Az. 
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To compare theory and experiment, we have adjusted 
the gain parameter A in order to get the best agreement. 
As mentioned on the curves of Fig. 6, we have shghtly 
modified the gain parameter replacing A by 1.1 x A for 
Az = 4 and 2 cm, and by 1.25 x A for Az = 10 cm. This 
small change of the gain parameter is interpreted here as 
small drifts (laser power,...) yielding small changes of the 
experimental signal and/or detection efficiency during 
the several hours of data acquisition. We have displayed 
the theoretical and experimental profiles on Fig. 6. The 
agreement is excellent. 

The Spd{z) curves displayed on Fig. 4 and Fig. 6 ex- 
hibit two regimes. The first regime is observed when 
Az is smaller than the width of the cloud of photons 
of Fig. 5 i.e. for Az < 1.5 cm. The width of the curve 
remains nearly constant, while the curve maximum is al- 
most proportional to Az. This is expected since in that 
case the width of the cloud is larger than PSF (i.e. for 
Az > 1.5 cm). The shape of the curves corresponds then 
to the shape of the cloud along z. On the other hand, the 
AOCT signal (i.e. the area of the Fig. 3 curves) is ever 
proportional to the number of diffuser that contribute 
coherently to the AOCT signal (i.e. is proportional to 
Az). As a consequence, since the width remains almost 
constant (and is equal to the width of the cloud), the 
maximum is proportional to Az. 

The second regime is observed when Az is larger than 
the width of the cloud, i.e. for Az > 1.5cm. The curve 
maximum saturates, and the width increases linearly 
with Az. This is expected since in that case the width 
of the cloud is smaller than PSF (i.e. for Az > 1.5 cm). 
One can then neglect the width of the cloud, and the 
Fig. 3 curves corresponds to the width of PSF(2;), which 
increases like Az, and thus like the AOCT signal. The 
maximum of the signal then saturates. 

4. AOCT images along the x and z directions. 

AOCT has been used to obtain 2D images of different 
samples. The random phase jump time is kept con- 
stant. The scan along z is made by varying the AOCT 
delay r with z = rCus- The scan along x is performed 
by moving the location x of the PZT transducer with a 
motor. In order to reduce the noise, the 2D acquired data 
are low pass filtered in the Fourier space with a cutoff 
frequency kc. 

A. Image of the cloud of Photon: sample without ab- 
sorbing inclusion. 

To validate the analysis of the PSF made in section 3, we 
have consider first a sample without inclusion in order 
to image the cloud of photons itself. 

We have first imaged a sample without absorbing in- 
clusion. The experiment is made with a /i^ = 6 cm~^ 
scattering sample, whose thickness is Ly = 3.2 cm, and 
whose dimensions are Lx x Lz = 6.5 x 6.5 cm^. The op- 
tical power is Ps = 630 mW/1 cm^ on the sample, and 
Pr = 650 mW/0.32 cm^ on the GaAs crystal. Acous- 
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Fig. 7. AOCT image (x, z) of a scattering sample without 
absorbing inclusion (e = 3.2 cm, ijl'^ = 6 cm~^). 
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Fig. 8. Profiles of the AOCT photodiode lock-in signal 
of the sample of Fig. 7 along x (a) and z (b) directions. 



tic pressure is Pmax = ^.6MPa., and lock-in integration 
time is Tc = 100 ms. Steps along x and z are 5x = 0.65 
mm, and 5z = 0.65 mm yielding 100 x 100 pixels images. 
The jump time is T<^ = 1.9 /is yielding Az = 2.9 mm. 
The data are low pass filtered in the Fourier space. The 
spatial cutoff frequency is /Cc = 1 mm~^ for kx and for 
kz- 

Figure 7 shows the xz image of the sample. Since there 
is no absorbing inclusion, we get an image of the cloud of 
photons that travel through the sample and that reach 
the photodiode. The US beam 

es the cloud, and the image corresponds to a cut of the 
cloud along the US beam z axis, and the PZT motion 
X axis. As shown on Fig. 7 the cloud is symmetric by 
rotation around the y axis, exhibiting the same extension 
within the x and z direction as shown on Fig. 8. This is 
expected since the x and z extension (17.5 and 16.8 mm) 
is much larger than the resolution along x (FWHM = 1.2 
mm for the focused US beam) and z {Az = 2.9 mm). 

B. Quantitative study of the z resolution in the high 
resolution case: images and cuts made with a low 
scattering sample (ji'^ 0) with two absorbing in- 
clusions. 

To perform a quantitative study of the z resolution in 
the high resolution case {Az down to 1.4 mm), we con- 
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Fig. 9. Photographs of the non scattering sample (e = 
3.0 cm, ji'^ = 0) having 2 absorbing inclusions (diameter: 
3 mm, spacing: 2 mm). 
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Fig. 10. AOCT xz images of the sample of Fig. 9 with 
Az = 11.4 mm (a), 5.7 mm (b), 2.9 mm (c) and 1.4mm 



sidered a sample with two absorbing inclusions of known 
geometry. These inclusions were made by Agar Agar with 
the addition of India ink drops. They have a cylindrical 
shape (diameter: 3 mm, length 5 mm), separated by 2 
mm, and oriented along y direction (see Fig. 9). In prat- 
ice, the absorption coefficient was not properly measured 
but was considered as high. We considered a low scat- 
tering sample with Agar Agar and without intralipid. 
This has been done to get a better sensitivity at high 
resolution, (since the AOCT signal is proportional to 
Az) an to we avoid shadow effect. Indeed, the density 
of photons near an absorbing inclusion is lower than ex- 
pected, because the diffused photons that travel through 
the sample may have crossed or will cross the absorbing 
zone. The scattering coefficient is nevertheless non 
zero (as seen on the image of Fig. 9 ) yielding a pho- 
ton diffusion cloud, whose size is not well known. Since 
AOCT is based on the detection of diffuse light, a dif- 
fuser (sheet of white paper) is placed at the outlet of 
sample. The experimental conditions are: Ps = 1 W/ 
1 cm^ on the sample; Pr = 630 mW/0.32 cm^ on the 
crystal; Pmax = 1.6 MPa; Tc = 100 ms; kc = I mm~^. 

We have studied the z resolution by imaging the sam- 
ple and by plotting vertical cuts for Az varying from 
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Fig. 11. Profiles of the AOCT photodiode lock-in signal 
(crosses and blues curves) of the sample of Fig. 9 along 
z direction with Az = 11.4 mm (a), 5.7 mm (b), 2.9 
mm (c) and 1.4mm (d), and theoretical profiles (dashed 
line curves) calculated with totaly (black ) and partially 
(red) absorbing inclusions. 



11.4 mm to 1.4 mm. The images are displayed on Fig. 10 
with Az = 11 A mm (a), 5.7 mm (b), 2.9 mm (c) and 
1.4 mm (d). Figure 11 presents their axial profiles along 
the vertical axis connecting the two inclusions (crosses 
and solid blue curves). On the image corresponding to 
Az =11.4 mm (a), we cannot distinguish any inclusion. 
We just see a decrease of the signal along the vertical 
column. For Az = 5.7 mm (b), the inclusions are seen, 
but the contrast is low. Changing to Az = 2.9 mm (c), 
the two inclusions are clearly resolved, but the contrast 
is not yet optimum. The latter is maximum for Az = 
1.4 mm (d), but the signal becomes quite low, especially 
near the absorbing inclusion. Image (d) is thus noisy. 
Typically, the best compromise between resolution and 
amplitude is obtained when the resolution is compara- 
ble to the dimensions of the object, that is to say for 
Deltaz = 2.9mm compared to the diameter of the inclu- 
sion of 3mm. We can also notice from image (a) to (d) 
that the acousto-optic signal decreases linearly with Az. 

In order to compare the experimental results to the- 
ory, we have also plotted on Figure 11 the theoreti- 
cal profiles Spd{z) (dashed curves). The calculation of 
Spd{z) = gl{z) (g) Cz{z) is made by considering that the 
density of photons Cz[z) within the inclusions is either 
zero (black dashed curves) or attenuated by a factor 0.3 
(red dashed curves) with respect to the density within 
the cloud. Moreover, the gaussian profile of the cloud 
(exp(— x^/k;^)) is supposed to be a little narrower than 
in section 3 (wg = 8.0 in place of 10.2 mm). As seen on 
Fig. 11, the imaging contrast obtained in experiment is 
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Fig. 12. Pictures of the diffusing sample (e = 2 cm, 
/i^ = 10cm~^) with two absorbing (3 mm diameter, 2 
mm spacing) at mid-thickness in the plane {yz). 
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Fig. 13. AOCT xz images of the sample of Fig. 12 with 
Az = 11.4 mm (a) , 5.7 mm (b) , 2.9 mm (c) and 1.4 
mm (d). 



lower than the one expected for totally absorbing inclu- 
sions (black dashed curves). This effect can be related to 
the imperfect absorption of the inclusion. By considering 
a finite absorption of 0.3, the agrement of the experiment 
with theory is much better (red dashed curves). 

C. Sample with two absorbing inclusions 

To illustrate how Az affects the image quality on a more 
realistic situation, we have imaged a scattering sample 
with two absorbing inclusions, whose are similar to those 
used in section B. These inclusions were made by the 
same mix of Agar Agar+Intralipides than the rest of the 
sample, with the addition of India ink drops. They have 
a cylindrical shape (diameter: 3 mm, length 5 mm), sepa- 
rated by 2 mm, and oriented along y direction (see Fig. 9). 
The sample thickness is L^y = 2 cm, with a scattering co- 
efficient /ig = 10cm~^. Two inclusions of 3mm diameter 
are oriented along the y direction and are spaced 2 mm 
(see Fig. 12). The experimental conditions are: Ps = 350 
mW/ 1 cm^ on the sample; Pr = 650 mW/0.32 cm^ 
on the crystal; Pmax = 1-6 Mpa; Tc = 100 ms; kc = I 
mm~^. 




20 < -20 -10 10 



z (mm) 
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Fig. 14. Profiles of the acousto optic AOCT photodiode 
lock-in signal of the sample of Fig. 12 along z direction 
with Az of 11.4 mm (a), 5.7 mm (b), 2.9 mm (c) and 1.4 
mm (d). 



Figure 13 shows the experimental images obtained for 
Az = 11.4 mm, 5.7 mm, 2.9 mm and 1.4 mm. Figure 
14 shows the corresponding vertical cut along z^ at the 
center of the sample. For Az =11.4 mm, we do not dis- 
cern any inclusion, but only a decrease in the signal at 
the center of the sample that informs us that it is prob- 
ably one or more absorbing inclusion. For Az = 5.7mm, 
both inclusions are discernible. This is accentuated for 
Az = 2.9 mm and 1.4 mm. We can briefly compare these 
images in scattering medium to those performed in the 
non-diffusing in previous section B. 

For Az = 2.9 mm, the inclusions are less visible here 
with scattering (/i^^ = 10cm~^), than without scattering 
(/i^ = 0). This can be interpreted as a shadow effect that 
can occur in scattering medium. Here again, we can also 
notice that the signal decreases linearly with Az. 

As in the section B, the best compromise between res- 
olution and amplitude is obtained when the resolution 
is comparable to the dimensions of the object, that is to 
say for Deltaz = 2.9mm compared to the diameter of 
the inclusion of 3mm. 

5. Conclusion 

In this paper, we have confirmed that the Lesaffre et 
al. [23] AOCT technique is able to bring z resolution to 
acousto optic imaging. We have presented new experi- 
mental results, and we have studied how the z resolu- 
tion depends on the jumping time T^, or on the corre- 
lation length Az. By working at low resolution, we have 
performed a quantitative comparison of the measured z 
shape with the theoretical one, and we have shown by 
this way that the theoretical expression of the AOCT 
Point Spread Function of reference [24] agrees with ex- 
periments. 
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We have imaged samples with 2 absorbing inclusion 
(3 mm diameter, 2 mm spacing) embedded within 2 dif- 
ferent samples (diffusing or not). The non diffusing sam- 
ple has been used to analyse the resolution and to com- 
pare it quantitatively with theory in the high resolution 
regime {Az = lAand2.9mm). The experimental resolu- 
tion agrees with theory. The scattering sample has been 
considered to get qualitative images closer to what ex- 
pected in living sample. The resolution varies similarly 
with Az but is slightly lower. We have interpreted this 
decreased resolution as a shadow effect. 

The non-scattering sample permits to validate the 
model by avoiding the problem of shadow effect which 
will tend to decrease the resolution. The corresponding 
images permit a quantitative analysis of the resolution in 
the case of a weak scattering sample with a good signal 
to noise ratio. Otherwise, the study with a high scat- 
tering sample show more realistic results, but in return 
more qualitative because the resolution is limited by the 
signal to noise ratio and the shadow effect. 

Note that the experiments have been done here with 
Agar -h Intralipid phantoms, and not with ex- vivo tumor 
tissues because phantoms are easier to prepare, and be- 
cause their optical properties (absorption and reduced 
scattering coefficients fi'^) are more reproducible, and 
can be thus more precisely measured. The experimen- 
tal analysis done here correspond to AOCT experiments 
made with photorefractive detection of the acousto op- 
tic signal [23] , which have been theoretically analyzed by 
Lesaffre et al. [24]. We did not try to analyse in detail 
the recent AOCT experiments made with holographic 
detection of the acousto optic signal [27] because the 
corresponding theory has not been published yet. More- 
over, the comparison of the figure of merit of the two 
detection schemes (photorefractive versus holographic) 
in a clinical in vivo application has not been done. 

In order to use the AOCT technique for biomedi- 
cal applications, as for example to get information on 
breast tumors, it is necessary to work with thicker scat- 
tering samples (> 5 cm), while respecting the acoustic 
biomedical standards. In practice, this is very difficult 
to realize at A = 1064 nm because of water absorption. 
On the other hand, experiments realized within the op- 
tical therapeutic window - i.e. around A = 800 nm - 
are easier because the absorption is lower and therefore 
the transmitted signal much bigger. This change of op- 
tical wavelength involves the use of other photorefrac- 
tive crystals (as SPS:Te, ZnTe, BSO) whose develop- 
ment is in progress [21,28]. Besides, the biomedical safety 
recommendations concern both the optical and acous- 
tical powers [29,30]. For the optical power, the limit 
is IW per cm at A = 1064nm for optical continuous 
emission, and all the experiments in this paper follow 
the standard. The acoustic power standards concern the 
peak to peak acoustic power and the average acoustic 
power. The first one requires a maximum instantaneous 
power of 190 VK/cm^, corresponding to a maximum pres- 
sure of 2.5MPa in water. All the experiments in this 



paper follow this first acoustic standard. The second 
acoustic standard requires an average acoustic power of 
720m W/cm^. Because of the acoustic continuous emis- 
sion of this set-up, this second acoustic is not reached. 
To satisfy this last standard, one could use the AOCT 
technique with long acoustic pulses (about 1 ms), which 
are compatible with the long pulse lasers already used in 
Ultrasound- modulated optical imaging [22]. These two 
modifications would allow the AOCT technique to be 
used with thicker samples, while respecting the acoustic 
biomedical standards. 

This work was supported by funds from the French 
National Research Agency (ANR-2011-BS04-017-SIMI 
4-ICLM). 
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